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Abstract. We study a nonlocal diffusion operator in a bounded smooth domain pre- 
scribing the flux through the boundary. This problem may be seen as a generalization of 
the usual Neumann problem for the heat equation. First, we prove existence, uniqueness 
and a comparison principle. Next, we study the behavior of solutions for some prescribed 
boundary data including blowing up ones. Finally, we look at a nonlinear flux boundary 
condition. 



1. Introduction 

The purpose of this article is to address the Neumann boundary value problem for a 
nonlocal diffusion equation. 

Let J : ~R N —>■ R be a nonnegative, symmetric J(z) = J{—z) with j RN J(z) dz = 1. 
Assume also that J is strictly positive in B(0, d) and vanishes in R^ \B(0,d). Equations 
of the form 

(1.1) u t (x,t) = (J * u — u)(x,t) = / J(x — y)u(y,t)dy — u(x,t), 

and variations of it, have been recently widely used to model diffusion processes, see pQ, 
0, P"Tj . [T3]. As stated in [7j if u(x, t) is thought of as a density at the point x at 
time t and J(x — y) is thought of as the probability distribution of jumping from location y 
to location x, then J RN J(y — x)u(y, t) dy = (J*u)(x, t) is the rate at which individuals are 
arriving at position x from all other places and — u(x, t) = — J RN J(y — x)u(x, t) dy is the 
rate at which they are leaving location x to travel to all other sites. This consideration, in 
the absence of external or internal sources, leads immediately to the fact that the density 
u satisfies equation (jl.lj) . 

Equation so called nonlocal diffusion equation, shares many properties with the 

classical heat equation w t = Am such as: bounded stationary solutions are constant, 
a maximum principle holds for both of them and, even if J is compactly supported, 
perturbations propagate with infinite speed. 

Given a bounded, connected and smooth domain, Q, one of the most common boundary 
conditions that has been imposed to the heat equation in the literature is the Neumann 
boundary condition, du/drj(x,t) = g(x,t), x £ dQ. 
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for x G f2. In this model we have that the first integral takes into account the diffusion 
inside Q. In fact, as we have explained the integral J J(x — y)(u(y,t) — u(x,t)) dy takes 
into account the individuals arriving or leaving position x from other places. Since we are 
integrating in Q, we are imposing that diffusion takes place only in Q. The last term takes 
into account the prescribed flux (given by the data g(x,t)) of individuals from outside 
(that is individuals that enter or leave the domain according to the sign of g). This is 
what is called Neumann boundary conditions. 

Our first result for this problem is the existence and uniqueness of solutions and a 
comparison principle. 

Theorem 1.1. For every uq G L x (Q) and g G L^ c ((0, oo); L 1 (R JV \ Q)) there exists a 
unique solution u of (jl.2|) such that u G C([0, oo); L l {VL)) and u(x, 0) = Uq(x). 
Moreover the solutions satisfy the following comparison property: 



Once existence and uniqueness of solutions is proved an important aspect in evolution 
equations is the asymptotic behavior as time evolves. In this context, we study the 
asymptotic behavior of solutions for certain fluxes on the boundary. 

First, we deal with a flux independent of time, that is, g(x,t) = h(x). As happens for 
the heat equation, in this problem, when h verifies a compatibility condition, we prove 
that solutions converge exponentially fast as t — > oo to the unique stationary solution of 
the problem with the same total mass as uq. If the compatibility condition is violated 
then solutions become unbounded as t — > oo. We have the following result. 



(1.3) 




Theorem 1.2. Let in addition J G L 2 (R N ). Let h G L 1 (R JV \ Q) such that 




Then there exists a unique solution <p of the problem 




that verifies J n u = f n ip and the asymptotic behavior of solutions of (jl.2|) is described 
as follows: there exists (3 = /3(J,Q) > such that 

(1.6) IK*) - ^IU 2 (^) < e -/3 *||w - <p\\L*(n)- 



If (jl-4j) does not hold then solutions of p. 2)1 are unbounded. 
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Next, we prescribe the boundary flux in such a way that it blows up in finite time. We 
consider a flux of the form 

(1.7) g(x,t) = h(x)(T-t)- a , 

with a nonnegative and nontrivial function h. 

For this problem we analyze the possibility that the solution becomes unbounded at 
time T a phenomenon that is known as blow-up in the literature. For blowing-up solutions 
we also analyze the rate of blow-up (that is the speed at which solutions go to infinity at 
time T) and the blow-up set (that is the spatial location of the singularities). 

We find that blow-up takes place in strips of width d (recall that J is positive in -5(0, d) 
and zero outside) around the support of h with blow-up rates that increase as the strips 
get closer to the support of h. 

Before stating our theorem we need some notation. We set Qo = Q, Bo = supp(/t) and 
define recursively for i > 1 

Bi = {x e £1 \ Uj<jSj : d(x, Bi-i) < d} 

and 

&i = &>i-l \ Bi. 

We also define the functions w^, Wi : R — > K by 



w i( x ) = 7 7T / J ( x - y)h(y) dy, 

[a - 1 



Wi(x) = - — — - / J(x - y)w i - 1 (y) dy for 1 < i < a, 
(a - 1) J r n\ Ui 

wi(x) = J(x- y)h(y)dy 



and 

Wi(x) = (a — l)tUi(x) for 1 < i < a. 

We can now state our result. 

Theorem 1.3. Let m addition J e L°°(R N ). Assume h G L°°(R N h > 0, 

In Ir n \u J( x ~ y)h(y) dydx ^ 0. Then, the solution of (jl.2|) with g(x, t) = h(x)(T — t)~ a 
blows up at time T if and only if a > 1 . 
If a > 1 is not an integer the blow-up set, B{u), is given by 



B(u)= |J B, 



Ki<\a] 



with the asymptotic behavior 

(T — t) a ~ l u(x, t) — > Wi(x) uniformly in Bi as t — > T 
for each i such that 1 < % < [a\. 
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If a is an integer the blow-up set, B(u), is given by 

B(u)= |J B u 



Ki<a 



with the asymptotic behavior, 

(T — t) a ~ l u(x,t) — ► Wi(x) uniformly in Bi as t 
for each i such that 1 < i < a and 

u(x, t) 



- ln(T - t) 



w a (x) uniformly in B a as t — ► T. 



Observe that blow-up in the whole domain (global blow-up) is possible. Indeed this 
happens for large values of a (depending on Q, h and d). 

One can compare this result with the corresponding one for the heat equation with 
boundary flux du/drj(x,t) = h(x)(T — t) a . For the heat equation solutions blow up if and 
only if a > 1/2 and in this case max^ u(x, t) ~ (T — £)~ Q+1 / 2 . Therefore the occurrence of 
blow-up and the blow-up rate for non-local diffusion are different from the corresponding 
ones for the heat equation. 

Finally we consider a nonlinear boundary condition of the form 

(1.8) g(y,t) = v?(y,t) 

where u is the extension of u from the boundary to the exterior of the domain in the 
following form: let us assume that a neighborhood of width d of dfl in M> N \ Q can be 
described by coordinates (z, s) where z G dQ and s is the distance from the point to 
the boundary, then we set u(z,s) = u{z). For this nonlinear boundary condition with 
nonlocal diffusion we have the following result. 

Theorem 1.4. Let in addition J G C(M ). Then, positive solutions blow up in finite 
time if and only if p > 1. As for the blow up rate, there exist constants c > and C such 
that 

(1.9) c(T - t)- 1/{p - 1] < maxu(M) < C{T - tJ-Vfr-i). 

X 

Moreover, the blow-up set is contained in a neighborhood of dVt of width K d, where 
K=[p/{p-l)\. 

There is a large amount of literature dealing with blow-up for parabolic equations 
and systems see for example the survey jU], the book ^H] and references therein. When 
blow-up is due to nonlinear boundary conditions see for example [TU], [12] , the surveys 
[3], jH] and the references therein. It is known that solutions of the heat equation with 
a nonlinear boundary condition given by a power blow up in finite time if and only if 
p > 1, the blow-up rate is given by ||w(a;,t)||i,«>(fi) ~ (T — t) _1 /( 2 (p -1 )) and the blow-up 
set is contained in dQ. Hence the blow-up rate and set are different but the blow-up set 
contracts to the boundary as the support of J becomes smaller. Observe that for J fixed 
the blow up set can be the whole domain Q if p is sufficiently close to 1 . 
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Organization of the paper. In Section El we prove existence, uniqueness and the 
comparison principle, in Section El we deal with the problem with g(x,t) = h(x), in 
Section 13] we analyze the blow-up problem and finally in Section El we study the problem 
with a nonlinear boundary condition. 



2. Existence, uniqueness and a comparison principle. 

In this section we prove Theorem 11.11 and give as remarks several consequences of the 
proof that will be used later in the paper. 

As in [5J, existence and uniqueness will be a consequence of Banach's fixed point theorem 
so we give first some preliminaries. 

Fix to > and consider the Banach space 



x to = c([o,t };L\n)) 



with the norm 



w 



max \\w{-,t)\\ L i m . 

(J<r<ro 



We will obtain the solution as a fixed point of the operator T : X to — » X tQ defined by 



(2.1) 



T W0:9 (w)(x,t) = w (x) + / J(x-y) (w(y, s) — w(x, s)) dyds 

Jo Jn 

/ J (x -y)g(y,s)dyds. 
Iq Ju N \n 



The following lemma is the main ingredient in the proof of existence. 

Lemma 2.1. Let w , z E L X (Q) ; g, h e L°°((0, * ); L X (R N \ Q)) and w,z E X tQ , then 
there exists a constant C depending only on Q and J such that 



IT 



u>0,9 



(w) - T ZQth (z)\\\ < \\w - zo\\ip-(q) + Cto||||w - 2||| + \\g - MU-tfo.to)^ 1 ^^))}- 



Proof. We have 



T Wo , g (w)(x,t) - T Zo , h (z)(x,t)\ dx < / | w - z Q \(x) dx 
n Jn 
t 



+ 



+ 



n 



o Jn 

t 



J{x-y) (w{y, s) - z(y, s)) - (w(x, s) - z(x, s)) 



dy ds 



dx 



o Jm. N \n 



J(x-y) {g(y, s) - h(y, s)) dyds 



dx. 
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Hence 

\T W0>g {w)(x,t) - T Z0:h (z)(x,t)\ dx < \\w - z \\ L i (n) 

+ / / \(w(y,s) - z(y,s))\dy+ / / \(w(x, s) - z(x, s))\ dx 
Jo Jn Jo Jn 

+ \g(y,s) - h{y,s)\ dydsdx. 

Jo Jw. N \n 

Therefore, we obtain, 

\\\ T w ,g( w ) ~ T zo,h( z )\\\ < Iko - z \\ L i { n) + Ct Q {\\\w - z\\\ + \\g - h\\ Lac{(0:to) . L i {R N\ n)) y 
as we wanted to prove. □ 

Theorem 2.1. For every uq € L X (Q) there exists a unique solution u of ()1.2j) such that 
u G C([0, oo); Moreover, the total mass in Q verifies, 

(2.2) / u(y,t)dy= / u (y) dy + / / / J(x - y)g(y,s) dydxds. 



o Jn 



Proof. We check first that T U0}9 maps X to into X t(i . From (j2.1|) we see that for < t\ < 

\\Tu ,g{w){t 2 ) - T Uihg (w)(ti)\\ L i { n) < 2 / / \w(y, s)\ dxdyds+ 

Jt x Jn 

+ / / \g(y,s)\dxdyds. 
Jh Ju N \n 

On the other hand, again from (J2.ll) 

\\T U0 ,g{w)(t) - W \\ L i {n) < Ct{||M|| + ||5 , ||L-((0,i );L 1 (M^\n))}- 

These two estimates give that T U0i9 (w) G C([0, t ]; L 1 ^)). Hence T UQt9 maps X to into X to . 

Choose t such that Ct < 1. Now taking zq = Wq = Uq, g = h in Lemma 12.11 we 
get that T uo g is a strict contraction in X <0 and the existence and uniqueness part of the 
theorem follows from Banach's fixed point theorem in the interval [0, to]- To extend the 
solution to [0, oo) we may take as initial data u(x,t ) G L 1 ^) and obtain a solution up 
[0, 2to\. Iterating this procedure we get a solution defined in [0, oo). 

We finally prove that if u is the solution, then the integral in Vt of u satisfies (|2.2jh 
Since t 

u(x, t) — uq(x) = / / J (x — y) (u(y, s) — u(x, s)) dy ds 



Jo Jn 

+ 11 J (x -y)g(y,s)dyds. 
Jo Jm N \n 

We can integrate in x and apply Fubini's theorem to obtain 

u(x,t)dx — / u (x)dx = I J {x — y) g{y, s) dy dx ds. 

Jn Jo Jn 
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and the theorem is proved. □ 

Now we give some consequences that we state as remarks for the sake of future refer- 
ences. 

Remark 2.1. Solutions of ()1.2)) depend continuously on the initial condition and bound- 
ary data. Let u be a solution of (J1.2)) with initial datum u and v a solution of (jl.2|) 
with g replaced by h and initial datum v . Then for every t > there exists a constant 
C = C(t ) such that 

max \\u(-,t) - v(-,t)\\ L i [n) < C\\u(-,0) - v(-,0)\\ L i (n) + C\\g - h\\ L oo mo) . L i {K N\ Q)) . 
Remark 2.2. The function u is a solution of (jl.2j) if and only if 

u{x,t)= e- Aix)t u {x)+ [ [ e- A{x){t ~ s) J(x-y)u{y,s)dyds 

(2.3) t J ° Jq 

+ [ [ e- A ^ t - s \j(x-y)g(y,s)dyds } 
Jo Jr n \q. 

where A(x) = f n J(x — y) dy. 

Observe that A(x) >a>0 (x G Q) for a certain constant a. 

Remark 2.3. From the previous remark we get that if u G L°°(Q x (0,T)), uq G C k (Q) 
withO <k<oo,ge L°°(R N \£lx(0,T)) and J G W k '\R N ), thenu(-,t) G C k (Ux [0, T])). 

On the other hand, if J G L°°(R N ), u G L°°(Q) and g G L\R N \ O x (0,T)), there 
holds that u G L°°(Q x (0,T)). (See Corollaru \2. r\ for an explicit bound in the case of 
continuous solutions). 

We now define what we understand by sub and supersolutions. 

Definition 2.1. A function u G C([0, T); L l ((Q)) is a supersolution of (jl.2|) ifu(x,0) > 
u (x) and 

(2.4) u t (x,t)> / J(x-y)(u(y,t)-u(x,t))dy+ / J(x - y)g(y,t) dy. 



Subsolutions are defined analogously by reversing the inequalities. 

Lemma 2.2. Let uq G C(Q), uq > 0, and u G C(Q x [0, T]) a supersolution to p. 2)1 with 
g > 0. Then, u > 0. 

Proof. Assume that u(x,t) is negative somewhere. Let v(x,t) = u(x,t) + et with e so 
small such that v is still negative somewhere. Then, If we take (xo,to) a point where v 
attains its negative minimum, there holds that to > and 

v t {x ,t ) = u t {x ,t ) + e > / J(x -y)(u(y,t ) -u{x ,t ))dy 



J(x- y)(v(y,t ) - v(x ,t Q )) dy > 
which is a contradiction. Thus, u > 0. □ 
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Corollary 2.1. Let J G L°°(R ). Let uq and t> L 1 (Q) wift -u > Vq and g, h G 
^((O.T);!, 1 ^ \ J])) wit/i g- > /i. Let w 6e a solution of $£2b with u(x,0) = u and 
Neumann datum g and v be a solution of ()1.2|) with v(x, 0) = v and Neumann datum h. 
Then, u>v a.e. 

Proof. Let w = u — v. Then, w is a supersolution with initial datum uq — vq > and 
boundary datum g—h > 0. Using the continuity of solutions with respect to the initial and 
Neumann data and the fact that J G L°°(W N ), we may assume that u, v G C(Q x [0,T]). 
By Lemma \2. 21 we obtain that w = u — v > 0. So the corollary is proved. □ 

Corollary 2.2. Let u G C(Q x [0,T]) (resp. v) be a supersolution (resp. subsolution) of 
JE2J). Tften, u > v. 

Proof. It follows the lines of the proof of the previous corollary. □ 

Corollary 2.3. Let u be a continuous solution of (II. 2j) with u(x, 0) = uq and Neumann 
datum g G L°°((R N \ O) X (0, T)) . Then, 

(2.5) u(x, t) < sup m + / sup / J(x — y)g(y, s) dyds. 

n Jo n Jrn\ q 

Proof. Let 

u(t) = supw + / sup / J(x — y)g(y,s) dyds. 
n Jo n Jrn\ q 

Then, v is a continuous supersolution of (jl.2|) . By the previous corollary we get the 
estimate (l2~5J) . □ 

Corollary 2.4. If J E L°°(R N ) and u G C([0, T]; L l (Q)) is a solution of JE2J) wift 
u(x, 0) G L°°(fi) ; # G L°° (0, T; L 1 (R iV \ ft)) £/ien, (|Z3I) /wZds. 

Proof. Let u„ be the solution of (jl.2|) with u n (x, 0) = ^o(^) and Neumann datum g„, G 
L°°([R N \Q)x (0, T)) such that g n -> # in L 1 ((R Ar \ fi) X (0, T)) and < -> w in L x (lT) 
with ||uo||L°°(n) < ||«o|U°°(n)- 

The result follows from the application of Corollary 12.31 to the functions u n G C(Q x 
[0, T]) and taking limits as n — > oo. □ 



3. Asymptotic behavior for g(x, t) = h(x). 

In this section we study the asymptotic behavior, as t — > oo, of the solutions of prob- 
lem ()1.2|) in the case that the boundary data is time independent. So we will assume 
throughout this section that g{x,t) = h(x) and that J G L 2 (R N ). We start by analyzing 
the corresponding stationary problem so we consider the equation 



(3.1) = / J(x - y){ V {y) - cp{x)) dy + / J(x - y)h{y) dy. 

Jn Jm. N \n 
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Integrating in f2, it is clear that a necessary condition for the existence of a solution (p 
is that 



(3.2) 0= / / J(x-y)h(y)dydx. 

Jn Jm. N \n 

We will prove, by means of Fredholm's alternative, that condition (J3.2)) is sufficient for 
existence and that the solution is unique up to an additive constant. 
To do this we write (|3.1|) in the form 

(3.3) <p(x)-K(<p)(x) = b{x), 
where 



b(x)=a(x) J(x -y)h(y)dy } 

Ju N \n 

K(tp)(x) = a(x) / J(x -y)(f(y)dy 



n 



and 



We consider the measure 



n 



a ( x ) = [ I J{ x ~ y) dy 

dx 



dfi 



a{x) 

and its corresponding space L 2 ^ of square integrable functions with respect to this measure. 

We observe that, due to our assumptions on J, the operator K maps L 2 into L 2 ^ and 
as an operator K : L 2 ^ — » L 2 ^ is compact and self adjoint. 

We look now at the kernel of / — K in L 2 . We will show that this kernel consist only 
of constant functions. In fact, let ip G ker(J — K). Then (p satisfies 



ip(x) = a(x) / J(x -y)<p(y)dy. 
Jn 

In particular, since J G L 2 (R N ), (p is a continuous function. Set A = maxip(x) and 
consider the set 

A = {x G H | (p(x) = A}. 

The set A is clearly closed and non empty. We claim that it is also open in Q. Let Xq G A. 
We have then 

(p(x ) = a(x ) / J(x Q -y)(p(y)dy. 



n 



Since a(xo) = (J n J(xo — y) dy)~ l and (p(y) < <p(xo) this implies (p(y) = <p(xo) for all 
y G Q fl B(xo,d), and hence A is open as claimed. Consequently, as Q is connected, 
A = Q and ip is constant. 

According to Fredholm's alternative problem (|3.1|) has a solution if and only if 

b(x)^- = 
alx) 
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or equivalently 

/ / J(x — y)h(y) dy dx = 0. 
Jn Ju N \n 

We have proved 

Theorem 3.1. Problem (j3.1j) has a solution if and only if condition (13. 2 j) holds. Moreover 
any two solutions differ by an additive constant. 

We will address now the problem of the asymptotic behavior of the solution of (jl.2|) . 
The next proposition shows the existence of a Liapunov functional for solutions of (jl.2|) . 
Its proof is a direct computation and will be omitted. 



Proposition 3.1. Let u(x,t) be the solution of (II. 2|) . Let us define 

1 

4 



F(u)(t) = — / / J(x — y)(u(y,t) — u(x,t)) 2 dy dx 



(3.4) 

J(x — y)h(y)u(x, t) dydx. 




n JR N \n 



Then 

g-F(u)(t) = -2 I {u t f{x,t)dx. 

We are now in a position to state and prove a result on the asymptotic behavior of 
continuous solutions. 

Theorem 3.2. Letu be a continuous solution of (jl.2)l with g(x,t) = h(x) where h satisfies 
the compatibility condition ()3.2j) . Let ip be the unique solution of (|3.1jl such that 



(f(x) dx = / uo(x) dx. 
Jn Jn 

Then 

(3.5) u(x,t) — > ip(x) as t — > oo 

uniformly in Q. 

When ()3.2|) does not hold, solutions of (jl.2|) are unbounded. 

Proof. Set w(x, t) = u(x,t) — <f(x). Then w satisfies 

w t {x,t) = / J(x-y)w(y,t)dy-w(x,t) / J{x-y)dy 



and J n w(x, t)dx = 0. 

By the estimate given in Corollary 12.31 we have that ||w||i°o(nx[o,oo)) is bounded in 
Q x [0, oo) by ||« - ^Hl-^)- 

Setting A(x) = J Q J(x — y) dy and integrating, the above equation can be written as 



w(x,t) = e- A(x)t w(x,0)+ / e - A{x){t - s) / J{x-y)w{y,s)dyds. 



n 
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We note that A(x) is a smooth function and that there exists a > such that A(x) > a 
for all x E Q. We observe that for x\, 12 £ ^ one has 

\ e -A( Xl )t _ e -A( X2 )t\ < e -«t t |^( Xl ) _ A{x 2 )\. 

With this inequality in mind it is not difficult to obtain, via a triangle inequality argument, 
the estimate 

\w(xi,t) - w(x 2 ,t)\ < D (\A(xi) - A(x 2 )\ + \w(xi,0) - w(x 2 ,0)\) 

where the constant D is independent of t. This implies that the functions w(-,t) are 
equicontinuous. Since they are also bounded, they are precompact in the uniform conver- 
gence topology 

Let t n be a sequence such that t n — > 00 as n — > 00. Then the sequence w(-,t n ) has 
a subsequence, that we still denote by w(-,t n ), that converges uniformly as n — > 00 
to a continuous function ip. A standard argument, using the Liapunov functional of 
Proposition 13 .![ proves that ip is a solution of the corresponding stationary problem and 
hence ip is constant. As j n w(x,t)dx = this constant must be 0. Since this holds for 
every sequence t n , with t n — * 00, we have proved that t) — > uniformly as t — » cxd as 
we wanted to show. 

When ()3.2|) does not hold the equation satisfied by the total mass, (12. 2|) . implies that 
m is unbounded. □ 

We end this section with a proof of the exponential rate of convergence to steady states 
of solutions in L 2 . This proof does not use a Lyapunov argument. It is based on energy 
estimates. 

First, we prove a Lemma that can be viewed as a Poincare type inequality for our 
operator. 

Lemma 3.1. There exists a constant C > such that for every u £ L 2 (n) it holds 

I (u(x) — (u)) 2 dx < C / / J(x — y)(u(y) — u(x)) 2 dy dx, 
Jn Jn Jn 

where (u) is the mean value of u in Q, that is 

1 f 

(u) = —— / u{x) dx. 
I"l Jn 

Proof. We can assume that (u) = 0. Now let us take a partition of M. N in non-overlapping 
cubes, Tj, of diameter of length h. Using an approximation argument we can consider 
functions u that are constant on each of the cubes T,, u\x t = 0%. We will only consider 
cubes Ti such that Tj fl Q 7^ 0. For this type of functions we have to prove that there 
exists a constant C independent of the partition such that 

^2\Ti\a 2 < C^2^2 / J(x ~ y)dydx(ai - a k ) 2 . 
Recall that there exist a > and r > such that J(x — y) > o for any \x — y\ < 2r. 



12 C. CORTAZAR, M. ELGUETA, J.D. ROSSI, AND N. WOLANSKI 

If the centers of two cubes T, TL are at distance less than r and h < r we have 

\2 



/ / J(x - y)dydx(a,i - a k ) 2 > a\Ti\\T k \(ai 

JTi JT h 



Given Ti, Tk two cubes intersecting Q there exists a number I, depending only on Q and 
r but not on h, such that there exist a collection of at most i, not necessarily pairwise 
adjacent, cubes Tj x , . . . ,Tj t intersecting Q with Tj x = T, Tj e = T k and such that the 
distance between the centers of T Jm and Tj m+1 is less than r. Since all the involved cubes 
have the same measure, we have 

m||T fc |(a,-a fc ) 2 <2M^|T Jm ||T Jm+1 | 



(3.6) 



\m=l 
-1 



<— V/" f J(x — y) dydx (a jm - aj ) 2 . 



m=l J im J i m +i 

The intermediate cubes used in (j3.6|) corresponding to each pair T, T k can be chosen 
in such a way that no pair of cubes is used more that a fixed number of times (depending 
only on the diameter of Q and r) when varying the pairs T, T k . Therefore, there exists a 
constant C, depending only on J and Q but not on h, such that 



i fc i k 

On the other hand, as we are assuming that 

1^1^ = o, 



I/) dydx (a,i - a k ) 



we get 

E i^ii* - a *) 2 ^ 2 N E i T ^ a *) 2 

i k i 

and the result follows. □ 
Now let us take the best J— Poincare constant that is given by 

/ / J{x-y)(u(y)-u(x)) 2 dydx 

(3.7) 0= inf . 

u&L2{n) / (u(x)-(u)) 2 dx 

Jn 

Note that by Lemma l3~T1 3 is strictly positive and depends only on J and fl. 

Now let us prove the exponential convergence of u(x, t) to the mean value of the initial 
datum when the boundary datum vanishes, i.e., h = 0. 

Theorem 3.3. For every uq G L 2 (Q) the solution u(x,t) of ()1.2|) with h = 0, satisfies 
(3-8) \\u(-,t) - («o)||| 2( n) < e-^Wuo - («o)||| w 

ifere /? given by ()3.7|) . 
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Proof. Let 

H(t) = ^J(u(x,t)-(u )) 2 dx. 
Differentiating with respect to t and using (|3.7jl . recall that (u) = (m ), we obtain 

H'{t) = -- l l J(x -y)(u(y,t) - u(x,t)) 2 dydx < f (u(x, t) - (u )) 2 dx. 
* Jn Jn * Jn 

Hence 

H\t) < -pH(t). 

Therefore, integrating we obtain, 

H(t) < e-^F(O). 

As we wanted to prove. □ 



As a corollary we obtain exponential decay to the steady state for solutions of ()1.2j) 
with h 7^ 0. 

Corollary 3.1. For every uq G L 2 (Q) the solution of ([1.2)1 . u(x,t), verifies 
(3-9) \\u - <f\\ 2 L 2 m < e-^\\uo - <f\\ 2 L 2 iQ y 

Here <p is the unique stationary solution with the same mean value of the initial datum 
and P is given by (|3.7jl . 

Proof. It follows from Theorem 13.31 by considering that v — u — (p is a solution of ()1.2D 
with h = 0. □ 

4. Blow-up for g(y,t) = h(y)(T -t)~ a . 

Now we analyze the asymptotic behavior of solutions of ()1.2)1 when the flux at the 
boundary is given by 

g(y,t) = h(y)(T-t)- a 
with h > and f n j R N\ n J{x — y)h(y) dy > 0. We will also assume that the initial data, 
and hence the solution, is non negative. Throughout this section u(x,t) will denote the 
solution of (jl.2j) with boundary and initial data as described above. Also in this section 
we will assume, without loss of generality, that T < 1. This makes the quantity — ln(T— t) 
positive which helps to avoid overloading the notation. 

Throughout this section we will assume that J 6 L°°(M. N ) and we will use the notation 
introduced in the Introduction. 

First, we prove that a = 1 is the critical exponent to obtain blowing up solutions. 

Lemma 4.1. The solution u(x,t) blows up at time T if and only if a > 1. 
Proof. Set 

M(t) = / u(x,t)dx, 
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then one has 

M'it) = 1 / f J(x-y)h(y)dydx> ° 

{1 - t) a J n J R N\ n [1 ~ t) a 

Hence, if a > 1 M{t) is unbounded as t /* T and the same is true for the solution u(x, t). 
On the other hand, if a < 1 we consider the solution of the ordinary differential equation 

C 

= (T-t) a With Z ^ = Z °' 

that is a supersolution of our problem if C and zq are large enough. Since z(t) remains 
bounded up to time T, a comparison argument shows that so does u(x,t). □ 

Lemma 4.2. There exists a constant C such that for each integer i such that 1 < i < a., 
the solution u(x,t) verifies 

C 

u(x, t) < - . in if i^a 

and 

u(x, t) < — Cln(T — t) in if i = a. 

Proof. If a > 1 we have that 

z(t) ( 



(t - ty- 1 

is a supersolution to our problem for G\ large enough, therefore 

(4-1) U(x,t) < (T _^)a-l 111 

If a = 1 the argument can be easily modified to get 

u(x, t) < -Ci ln(T - t) in Q . 

Now for i G fli we have 
(4.2) ut(x,t)= / J(x -y)(u(y,t)-u(x,t))dy 

which implies 



(4.3) u t (x,t) < / J(x -y)(u(y,t) -u(x,t))dy+ J(x - y)u(y,t) dy 
Assume that a > 2. In this case, in view of (j4.1|) . we can use the function 

with C*2 large enough, as a supersolution in f2i to obtain that 

u(x, t) < — ^ — -, in fii. 
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As before if a = 2 we get 

u(x,t) < -C 2 ln(T-t), in Q v 

The previous argument can be repeated to obtain the conclusion of the lemma with the 

constant C = max C,-. □ 

i<i<[«] 

We can describe now precisely the blow up set and profile of a blowing up solution. 
Theorem 4.1. If a > 1 is not an integer the blow-up set, B(u), is given by 

B{u)= |J B t , 

l<i<[a] 

with the asymptotic behavior 

(T — t) a ~ l u(x, t) — > Wi(x) uniformly in Bi as t —>T 

for each i such that 1 < i < [a\. 

If a is an integer the blow-up set, B{u), is given by 

B(u)= |J B h 

l<i<a 

with the asymptotic behavior, 

(T — i) a ~ l u(x, t) — > Wi(x) uniformly in Bi as t — > T 

for each i such that 1 < i < a and 

u(x,t) fin i rp 

— ^ — — — > w a {x) uniformly in o a as t — > 1 . 

Proof. We have 

Kv) 



dy. 



(4.4) ut(x,t)= [ J(x-y)(u(y,t)-u(x,t))dy + J(x-y) 

Jn JR N \n K 1 ~ t) 

We prove first the theorem in the case when a — 1. Integrating ()4.4)1 in t and using 
that, by Lemma flOJ), u(x,t) < -Cln(T - t) we get 



u(x, t) 



ln(T - t) 
This proves that 



u(x,ti) „ 1 f . . InT 



lim — j ^ - ^ = uniformly in fi . 



Also if x e fii (|Ojl reads 

M 4 (x,t)= / J(x -y)(u(y,t) -u(x,t))dy. 
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Integrating in t and using again Lemma ()4.2j) we have 

u(x, t) < u(x, 0) — C / In (T — r)dr. 



o 



Hence u is bounded in VLi and the theorem is proved if a = 1. 
Assume now that a > 1 and consider the change of variables 

v x {x, s) = {T- t) a ~ l u{x, t), s = - ln(T - t). 

Since u verifies f!4.4|) . V\ satisfies 



(vi) s (x,s)= e s J(x-y)(v 1 (y, s) - v x (x, s))dy 
(4.5) Ja 

+ / J(x -y)h(y)dy- (a-l)vi(x,s). 



Integrating in s we obtain 

Vi(x, s) — Wi(x) = e~( a ~ lS)s v(x, 0) 

( 4 - 6 ) +e -(«-i)s f e {a ~ 2)r [ J(x-y)(v 1 (y,r)-v 1 (x,r))dydr 



Jo Jn 
- e -("-i)' Wl ( x ). 

If a 2 since, by the previous lemma, v i is bounded we get 

(4.7) Mx, s) - wi(x)\ < C{e~ s + e - {a ~ 1)s ) 
for some constant C. This implies that 

(4.8) (T -t) a ~ l u(x,t) -> 

uniformly in f2 as t — > T. 

We note that if a < 2, since w x (x) vanishes in f2 x , (|4.7j) implies 

K(x,s)| < C"e- (Q - 1)s for x e 

and hence 

u(x,t) <C for x G 

Consequently, if 1 < a < 2 the blow up set of u is Sl \ ^i = &i and the asymptotic 
behavior at the blow up time is given by ()4.8j) . 

We have to handle now the case a = 2 which is slightly different. In this case instead 
of estimate ()4.7|) there holds 

(4.9) \vi(x, s) - Wl (x)\ < C{se~ s + e~ {a - 1)s ). 
This still implies that 

(4.10) (T -t) a ~ l u(x,t) -> w x {x) 
uniformly in Qq as t — > T but does not ensure that u is bounded in VL\. 
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If x G Oi one has 

u t (x,t) = J(x-y)(u(y,t)-u(x,t))dy + J(x -y)u(y,t) dy - J(x -y)u(x,t) dy. 
Jn x JBx Jb 1 

Integrating in t we obtain 

u(x,t) - / J(x - y)u(y,r) dydr = u(x,0) + h + I 2 , 

JO JBi 

Ii(x,t) = / / J(x — y)(u(y,r) — u(x,r)) dy dr 
Jo Jili 



where 



and 



I 2 (x, t) = — / / J(x — y)u(x,r) dy dr. 
Jo JB^ 



'0 JBi 

Now using the fact that for z G Q\ one has u(z,t) < — Cln(T — t), it can be checked 
that 

j ^ ' ^ — > uniformly in as t — > T 

and also 

72 t) ^ ^ uniformly in fix as t — > T. 
ln(T — t) 

Moreover since (T — t) — * uniformly in B\ as t — > T one has 

- i /J" / / J{x-y)u(y,r)dydr^ [ J(x - y)w 1 (y) dy 

uniformly in fii as t — > T. 

Putting together this information we deduce that 

-\$^T)^ J Bi J (*-y)My)dy 

uniformly in Q± as t — > T. 

Finally since t) < — C ln(T — t) in Q± we can argue as in the proof of the case a = 1 
to show that u remains bounded in Q 2 . So we have shown that if a = 2, then the blow 
up set, B(u), of u is given by 

B{u) = B l UB 2 

with the asymptotic behavior 



and 

If a > 2 setting 



lim(T — t)u(x,t) = Wi(x) uniformly in Q 



lini(T — t)u(x, t) = w 2 (x) uniformly in Q 1 . 



v 2 (x, s) = (T - t) a ~ 2 u(x, t), s = - ln(T - t), 
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we obtain for x G Q\ the equation 



(4.11) (v 2 ) s (x,s) = e 

This can be written as 

(v 2 ) s (x, s) = e 

(4.12) 



J(x - y)(v 2 (y,s) - v 2 (x,s))dy - (a - 2)v 2 (x, s). 



J(x- y)(v 2 (y,s) - v 2 (x,s))dy 



+ / J(x-y)vx(y,s)dy-v 1 (x,s) / J(x-y)dy 

J Bi JBx 

-(a - 2)v 2 (x, s). 

Again integrating in s , after observing that by ()4.7|) \vi(x, s)\ < Ce~ s since x G we 
obtain that 

f s 



(4.13) 



v 2 {x,s)-e- {a ~ 2)s e - (a " 2)r / J{x-y)v 1 {y,r)dydr 



Bi 



< C(e~ s + e 



for some constant C provided that «^3. 
Also by ()4.7j) one has that 



J(x-y)vx(y,s)dy- J(x - y)w 1 (y) dy 



Bi 



< Ce 



and consequently 

w 2 {x) -e~ {a - 2)s [ e" (Q - 2)r f J(x-y) Vl (y,r)dydr <Ce- s . 

JO JBi 

This, together with f)4.13|) . implies that for all x G Qi 

\v 2 (x, s) - w 2 {x)\ < C{e~ s + e - {a ~ 2)s ) 
for some constant C and hence 

(T -t) a ~ 2 u{x,t) ^w 2 (x) 

uniformly in Qi. 

The above procedure can be iterated to obtain for all integers i such that 1 < i < a 

(4.14) \vi(x, s) - Wi (x)\ < C(e~ s + e - {a ~ i)s ) for all x G Q^i 
for some constant C and hence 

(4.15) (T - t) Q ^M(x, t) -> Wi(x) 

uniformly in tti-i. Moreover, it follows from ()4.14j) that for x G fl[ a ], if a is not an integer, 

\v [a] (x,s)\ <Ce-^~^ s 

and hence 

u(x,t) < C for all x G fl[ a }. 

In this fashion we have proved that, if a is not an integer, the blow up set of u is 
Ui<t<[a] &i and the behavior of u near time T in Bi is given by f)4.15j) . This proves the 
theorem in the case that a is not an integer. 
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In the case that a is an integer one can argue as in the proof of the case a = 2 to obtain 
the result in that case. □ 

5. Blow-up with a nonlinear boundary condition. 
In this section we deal with the problem 



u t = / J(x - y)(u(y,t) - u(x,t))dy + / J(x - y)u p (y,t)dy. 
(5.1) Jn Jm.N\ n 

u(x, 0) = Uq(x). 

Here we assume that J G C(R JV ), Uq G C(Q), u > and w is the extension of u to a 
neighborhood of Q defined as follows: take a small neighborhood of dfl in M. N \ Q in such 
a way that there exist coordinates (s, #) G (0, s ) x <9fi that describe that neighborhood 
in the form y = z + srj(z) where z G dfl and r](z) is the exterior unit normal vector to dQ 
at z. We set 

u(y,t) = u(z,t). 

We also assume that d < so therefore for any x G Q the ball centered at x and of radius 
d is contained in the above mentioned neighborhood. 

We address now the problem of local existence in time and uniqueness of solutions. 

As in the previous sections we set 

A(x) = / J(x- y)dy 
Jn 

and observe that there exists a > such that A(x) > a for all x G Vt. 

As earlier we obtain a solution of (j5-H) as a fixed point of the operator T defined by 

Tu(x,t)= e- A{x)t u {x) + [ e - A{x){t ~ s) [ J(x-y)u(y,s)dyds 

Jo Jn 

+ / e - A(x){t - s) [ J(x-y)u p (y,s)dyds. 
Jo Jvi N \n 

We split the proof of existence into two cases. We deal first with the case p > 1 since in 
this case we have uniqueness of solutions. In this direction we have the following theorem. 

Theorem 5.1. 

a) Let p > 1. There exists t Q > such that problem (15. lj) has a unique solution 
defined in [0, to). 

b) Let p < 1. There exists to > smc/i £/ia£ problem (j5.1|) /ias at /east one solution 
defined in [0, t )- 

Proof. Fix M > ||mo||oo, to > and set 

X = iu G C7(Hx [0,t )) /u > , = sup \u(x,t)\ <2M\. 

^ (x,t)enx[o,t ) 
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If to is chosen small enough, then T maps X into X. Indeed, we have for t < to and 

uex 

ft 



\Tu(x,t)\ < e- A{x)t u (x) + 



+ 



-A(x)(t-s) 



e -A( X ) { t-s) / j( X -y)\ U (y jS )\ dyds 

J{x -y)\vP{y,s)\dyds 



iR N \Q 

<M + t (2M + (2M)p) < 2M 

if to is small. 

Proof of a): We will prove that for p > 1 we can choose to in such a way that T is a 
strict contraction. In fact, for t < to and u± , u 2 G X 

\T Ul (x,t)-Tu 2 (x,t)\ < [ e- A ^ l - s) [ J(x-y)\u 1 (y,s)-u 2 (y,s)\dyds 



+ 



-A(x)(t~s) 



J(x-y)\u{(y,s) -u p 2 {y,s)\dyds 



<to(l+p(2My- 1 )\\\u 1 -u 2 \\\ 
and part a) of the theorem follows via Banach's fixed point theorem. 

Proof of b): We have that T maps X into X if to is small enough. We claim that the 
operator T : X — > X is compact. Indeed, for t\, t 2 < to, u G X and X\, x 2 G £1 we have 

\Tu{x u t x ) - Tu(x 2 ,t 2 )\ < e- A ^Uo{x x ) - e- A ^ t2 u (x 2 ) 



+ 



+ 



e -A( X1 )( tl -s) / J( Xl -y) U (y jS )dyd S 

Jn 

e -A{x 2 )(t 2 -s) 

Jn 

J(xi -y)u p (y,s) dyds 



J(x 2 - y)u(y, s) dyds 







-A(xi)(ti-s) 
-A(x 2 )(t 2 -s) 



J(x 2 -y)vP(y, s) dyds 



R N \Q 

As in the proof of Theorem 13.21 we have that for Xi, x 2 G Q one has 

| e -A(xi)t _ e -A0«2)t| < e -at f _ A ^ 

with a > 0. This inequality plus the fact that J is integrable imply, via a triangle 
inequality argument, that the family {Tu / u G X} is equicontinuous and, since it is 
bounded, it is precompact in (X, \ \\ ■ |||). Consequently, since T is clearly continuous in 
X, it is a compact operator and the claim is proved. Part b) of the theorem now follows 
from Schauder's fixed point theorem. □ 
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Remark 5.1. We observe that the same argument of the proof of part a) of Theorem \5.1\ 

provides existence of a unique solution if the boundary nonlinearity takes the form f(u) 
with f locally Lipschitz. 

Now we prove a comparison lemma for solutions of (|5.1|) . 

Lemma 5.1. Let u be a continuous subsolution and v be a continuous supersolution of 
problem (15. lj) defined in [0, to) ■ Assume u(x, 0) < v(x, 0) for all igO. Then 

u(x, t) < v(x, t) 

for all (x,t) G x [0, t ). 

Proof. Assume, for a contradiction, that the lemma is not true. Then, by continuity, there 
exists xi 6 fl and < ti < to such that u(xi,ti) = v(xi,ti) and u(x, t) < v(x,t) for all 
(x,t) G Q x [0,ti). We have now 

= u(x u ti) - vfa, h) = e ~ A{xi)tl (u(x u 0) - v(x u 0)) 

+ r e -^)(*i-) f J( Xl -y)( u (y,s)-v(y,s))dyds 
Jo Jn 

+ f 1 e -Mxi)(tis) f J( Xl -y)(u p (y,s) -v p (y lS ))dyds <0 
Jo Jr n \q 

a contradiction that proves the lemma. □ 

Now we use this comparison result to prove the lack of uniqueness for p < 1. 

Proposition 5.1. In the case p < 1 with uq = t/jere exists a nontrivial solution of 
problem fl5.1|) . Hence this problem does not have uniqueness. 

Proof. Let be a positive solution of b' = b p with 6(0) = and < a(x) < 7 be a 
continuous function with a(x) =7 on dQ. Let 7 > be so small as to have 



Y / J(x -y)dy> 2 7 
Ju N \n 

for every x G £1. Then, 

v(x, t) = a(x)b(t) 

is a subsolution to our problem for a certain interval of time, (0,t )- 

Let e > be given and consider a locally Lipschitz function f £ such that / e (s) = s p 

for s > e/2. It follows from Remark 15.11 that there exists a unique solution, w £ , of (j5.1j) 

with the boundary nonlinearity replaced by f £ {w) and initial data w £ (x,0) = e. By the 

comparison principle w £ > e and hence it is a solution of (|5.1|) . 

By comparison, the sequence w e is monotone increasing in e. In particular, for every e 

w £ is defined on the interval [0, ti] where W\ is. Therefore, by monotone convergence, we 

obtain that the limit 



w = lim w £ 

e^0 



is a solution with w(x, 0) = 0. 



22 



C. CORTAZAR, M. ELGUETA, J.D. ROSSI, AND N. WOLANSKI 



Using again the comparison principle we obtain that w £ (x,t) > v(x,t) for < t < 
mm{t ,ti}. Hence, w(x, t) > for every < t < min{t ,£i} and all x G supp(a). □ 

We address now the blow up problem for solutions of (|5.1|) . In this direction we have 
the following theorem. 

Theorem 5.2. 

a) Let p > 1, then every non trivial solution of (|5.1|) blows up in finite time. 

b) Let p < 1, then every solution of (|5.1j) is globally defined in time, by this we mean 
that it exists for all t G [0, oo). 

Proof. Proof of a): Let u be a solution of (|5.1|) and assume, for a contradiction, that it is 
globally defined in time. 

Since u > and L J(x — y) dy < 1 we have for x 6 Q, 



Ut(x, t) > —u(x, t) + 



J(x- y)vF(y,t) dy. 



Here we have used that the equation is satisfied for x G dfl. 

Integrating on dQ, denoting by dS x the surface area element of dQ, we get 

— / u(x,t)dS x > - / u(x, t) dS x + [ [ J(x - y)u p (y,t) dydS x . 
dt Jan J an Jan Jr n \u 

J(x — y)dydS x > 

J(x - y)u(y,t) dydS x 

an Jr n \u 



Since 



ion Jr n \u 

an application on Jensen's inequality implies that 

/ u(x,t) dS x > 
at Jan 

for some constant C > 0. 
Now, 



u(x, t) dS x + C 



J(x-y)u(y,t) dydS x 



an Jr n \d. 



/ / / J(x — a — srj(a) s ju(a,t) dS a ds dS x 
Jan Jo Jan 

j(x — a — sr](a)) ds u(a, t) dS a dS x 



an Jan 



>5\dQ\ [ u(a,t)dS 
Jan 



since we have that J(z) > c > for \z\ < d/2. 



Thus, if we call 



m(t) 



u(x,t) dS x 



an 
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we have 

(5.2) m'{t) > -m(t) + jm p (t). 

This implies that m(t) — > oo in finite time if for some t , m(t ) is large enough. 

Since we are assuming that u(x,t) is defined for every t > 0, it holds that m{t) is 
defined (and finite) for all t > 0. Let us see that this leads to a contradiction. Let v (x, t) 
be the solution of ()1.2|) with g = and v(x, 0) = u . By Theorem 13.21 we get that 



uniformly in Q. Since u is a supersolution for the problem satisfied by v, there exists 
ti > such that for t > t 1; 

> — — / M = c > 0. 
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Therefore, 

M(t) = / u(x,t)dx 




>M(t 1 )+ / / J(x -y)u p (y,s)dydxds 
'h Jn 



> M(t x ) + (t-t x )ccg. 
Arguing as before we get that 



«0M)> 2^T M (^) 



for i large enough. Now 



m(t ) = I u(x, t ) dx > ^M(t 2 ) > C(M{t x ) + (t 2 - h)c eg). 



<-)<> 



2IOI 



This implies that m(t ) is as large as we need if to is large enough, hence m(t) is not 
defined for all times and we conclude that u blows up in finite time. 

Proof of b): Let p < 1 and let u be a solution of (j5.1J) . Set t> (t) = C(t + l) 13 ^. It is 
directly checked that 



i 



v'(t) = VP (t . 

1 — p 



> max / J(x — ?/) c/y 



Picking C such that 

1 — V ~ xsTi 

we have that v is a supersolution of (|5.1j) . Moreover taking C larger, if necessary, such 
that m(x, 0) < f (0) in O we obtain by Lemma f5. II that 

u(x, t) < v(t) 
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as long as u is denned. This implies the theorem in the case p < 1. The case p = 1 is 
proved in the same fashion but using v(t) = Ce t as a supersolution. □ 

Our next result is an estimate of the blow up rate of blowing up solutions of (j5.1j) . 

Theorem 5.3. Let u be a solution of (|5.1|) that blows up at time T . Then there exists a 
constant C such that 

(5.3) (p_i)-V(p-i)( T _ f )-V(p-i) < \\ u (;t)\\L°°(n) <C(T-ty l ^~ l \ 
Proof. Let 

v (t) = { P - i)-v(p-i)(r - t)- 1/ip - 1) . 

One can easily check that v is a supersolution of our problem. If for some to G (0, T) one 
has 

IK->*o)IU°°(n) < u(to)> 
then there exists T > T such that 

<Mo) < (p- i)-V(p-i)(f -t )-v(p-i). 

Let v(t) := (p — 1) _1 /(p _1 )(T — i) -1 /^- 1 ) that is also a supersolution to our problem in 
the interval [t , T). 

Using a comparison argument, we obtain 

u(x,t) < (p-l)" 1 /^-!)^ _ t )-V(p-l) j 

for all £ G (t ,T). This contradicts the fact that T >T and hence 

(p _l)-i/( P -i) (T _ t )-i/( P -i) < ||<,*)|| ioo(n) . 

The proof of the reverse inequality is more involved. By equation (|5.2|) . if 
m(i) = / u(x,t) dS x — > oo as £ ^ T, 

we have 

(5.4) m(t) < C(T -t)- l l^- l \ 
Now, we claim that 



V(p-i) 



/ / J(x-y)v?(y,s)dyds<C, 
Jo Jm. N \n 



(5.5) (T - f) 

for all x G <9f2. In fact, if this does not hold, there exists a sequence (x n , t n ) with x n G dfl, 
t n /" T, such that 

(T - t„) 1/(p - 1} /"" / J(x„ - y)M p (i/, s) (fyds - oo. 
Jo JR N \n 

By compactness we may assume that x n — * Xq G <9f2. Hence 

(T-t n )V(p-*) /"" / u p (y,s)dyds -> oo. 

JO J(R N \£l)nB(x ,2d) 
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Therefore there exists a point xi G dfl such that for a subsequence that we still call t n , 

„ n t*n 

iV(p-i) 



(T — t n ) ,(p ' I / u?(y,s)dyds^oc. 

JO J (R N \Q)nB(x!,d/A) 

Since every function involved is nonnegative and J(z) > c > for \z\ < d/2 we get 

for every x G <9f2 Pi {\x — x%\ < d/4}. 
Using (|2.Hjl . we get 

( T _ tn )l/( P -l) M ( £)tn ) > c ( T _t n )V(p-l) /"" /" J( X -y)v?(y, S )dyds^OC. 

Therefore, 



JR^fi 



( T _ tn )i/(p-i) m ( tn ) = ( T _ tn )i/(^i) / u (x,t n )dS x 



which contradicts (|5.4jl . The claim is proved. 

Using again that J(z) > c > for 2 < d/2 we get that (|5.5jl holds for every x G fi. In 
fact, first we see that for every x G dVL 

(r-t)V(p-i) /" /" s ) dS a ds < C. 

Jo JdnnB d/4 (x) 

Then, since dQ is compact we deduce that 



( T _ t )i/(p-i) f 
Jo 



u p (o~, s) dS a ds < C. 
an 



This immediately implies, by using that J G L°°, that ()5.5|) holds for every 16H. 
Now, let for to < T, 

M = _max (T-ty/b-Vufat) = (T - ti) 1/&,_1) u(xi, 
nx[o,to] 

This implies by using again (|2.Hj) that 

M<C+ j 1 e ~M^){ti-s) f j( Xl - y )Mdyds <C + (l-e~ A{xi)tl )M. 
Jo Jn 
So that, since A(x) > a > 0, 

M < C, 

with C independent of to- The result follows. □ 

Corollary 5.1. Let u be a solution of (j5.1j) that blows up at time T. Then, the blow-up 
set, B{u), verifies 

(5.6) B(u) C {x G Q I dist (x, dQ) < Kd} 

where K = [p/ (p — 1)] . 
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Proof. The proof follows from the results in Section □ 
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